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ABSTRACT 

With close pairs of quasars at diflFerent redshifts, a background quasar sightline can be used to study a 
foreground quasar's environment in absorption. We search 149 moderate resolution background quasar 
spectra, from Gemini, Keck, the MMT, and the SDSS to survey Lyman Limit Systems (LLSs) and 
Damped Lya systems (DLAs) in the vicinity of 1.8 < z < 4.0 luminous foreground quasars. A sample 
of 27 new quasar-absorber pairs is uncovered with column densities, lO"'^^ '^ cm~^ < A^hi < lO^^'^cm"^, 
and transverse (proper) distances of 22 kpc < R < 1.7 h^^ Mpc, from the foreground quasars. If 
they emit isotropically, the implied ionizing photon fluxes are a factor of ~ 5 — 8000 times larger than 
the ambient extragalactic UV background over this range of distances. The observed probability of 
intercepting an absorber is very high for small separations: six out of eight projected sightlines with 
transverse separations R < 150 h^^ kpc have an absorber coincident with the foreground quasar, of 
which four have A'hi > lO^^cm"^. The covering factor of TVhi > lO^^cm"^ absorbers is thus ^ 50% 
(4/8) on these small scales, whereas < 2% would have been expected at random. There are many 
cosmological applications of these new sightlines: they provide laboratories for studying fluorescent 
Lya recombination radiation from LLSs, constrain the environments, emission geometry, and radiative 
histories of quasars, and shed light on the physical nature of LLSs and DLAs. 

Subject headings: quasars: general - intergalactic medium - quasars: absorption lines - cosmology: 
general - surveys: observations 



1. INTRODUCTION 

Although optically thick absorption line systems, that 
is the Lyman Limit Systems (LLSs) and damped Lyman- 
a systems (DLAs), are detected as the strongest absorp- 
tion lines in quasar spectra, the two types of objects, 
quasars and absorbers, play rather different roles in the 
evolution of structure in the Universe. The hard ultravi- 
olet radiation emitted by luminous quasars gives rise to 
the am bient extragalactic ultraviolet (UV) b ackground 
(see e.g. iHaardt fc Mad au 1996; Meiksin 2005) responsi- 
ble for maintaining the low neutral fraction of hydrogen 
(~ 10"^) in the intergalactic medium (IGM), established 
during reionization. However, high column density ab- 
sorbers represent the rare locations where the neutral 
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fractions are much larger. Gas clouds with column den- 
sities log A^Hi > 17.2 are optically thick to Lyman contin- 
uum (tll ^ 1) photons, giving rise to a neutral interior 
self-shielded from the extragalactic ionizing background. 
In particular, the damped Lya systems domina te the 
neutr al gas content of the Universe llProchaska et alJ 
|2QQ3)j which provides the primary reservoir for the star 
formation which occurred to form the stellar masses of 
galaxies in the local Universe. 

One might expect optically thick absorbers to keep a 
safe distance from luminous quasars. For a quasar at 
z = 2.5 with an r-band magnitude of r = 19, the flux 
of ionizing photons is 130 times higher than that of the 
extragalactic UV background at an angular separation of 
60", corresponding to a proper distance of 340 kpc 
and increasing as r^^ toward the quasar. Indeed, the 
decrease in the number of optically thin absorption 
lines (logTVni < 17.2 hence tll ^ 1), in t he vicinitv of 
quasa rs, known as the proximity effect ijBaitlik et alJ 
11988ft . has been detected and its st rength prov ides a 
measurement of the UV background (Sc ott et"al] l2000'l. 
If Nature provides a nearby background quasar sightline, 
one can also study the transverse proximity effect, which 
is the expected decrease in absorption in a background 
quasar's Lya forest, caused by the transverse ionizing 
flux of a foreground quasar. It is interesting that the 
transverse effect has yet to be detected, in spite of 
many atte mpts (iCrottsI [l989: PobrzYcki fc Bcchtold l 
1991 ; iFernandez-Soto. Barcons. Carballorfc Weba 
1 9951 ILiske fc Willigeil ^(iml 

Schir ber. Miralda-Escude. fc McDonald! 120041: iCrofli 



120041. but see Jakobsen et al.2003). 

On the other hand, it has long been 
known that quasars are associated with en- 
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hancements in the distrib u tion of galaxies 
jBahcall. Sch midt, fc ChlTil ITflfiflt IVee CreenI ITflSl 
1987 1 iBahcaU fc Chokshilll99lHSmith Bovle. fc Maddoxl 
200ntlBrown. Bovle. fc Websteill20niHSerber et al.ll2n0fil: 

Coil et al.l2 006'). although these measurements of quasar 
galaxy clustering are limited t o low redshifts < 1.0. 
Recently, Adclbcr ger fc Steidell l)2005j) . measured the 
clustering of Lyman Break Galaxies (LBGs) around lu- 
minous quasars in the redshift range (2 < z < 3.5), and 
found a best fit correlation length of tq = 4.7 h^-^ Mpc 
(7 = 1.6), very similar to t he auto-correl at ion length of 
z ^ 2 - 3 LBGs ( Adclbcrg er et al.ll2003j) . ICooke et all 
l)2006D recently measured the clustering of LBGs around 
DLAs and measured a best fit tq = 2.9 h^^ Mpc 
with 7 = 1.6, but with large unce rtainties (sec also 
IGawiser et all l200lt IBouche fc Lo^enth al 2004J. If 
LBGs are clustered around quasars, and LBGs are 
clustered around DLAs, might we expect optically 
thick absorbers to be clustered around quasars? This is 
especially plausible in light of recent evidence that DLAs 
arise from a high redshift gala xy population which are 
not unlike LBGs ijMoller et al.l l2002). 

Clues to the clustering of optically thick absorbers 
around quasars come from a subset of DLAs with Zabs 
Zcm known as proximate DLAs, which have absorber 
redshifts within 3000 k m s~^ of the emiss ion redshift 
of the quasars (see e.g. iMoUer et alJ I199E ). Re cently, 
iRussell et aP ()2005j) Csee also lEUison et al.ll200l . com- 
pared the number density of proximate DLAs per unit 
redshift to th e average num ber density of DLAs in the 
the Universe ijProchaska et al. 2005). They found that 
the abundance of DLAs is enhanced by a factor of ~ 2 
near quasars, which they attributed to the clustering of 
DLA-galaxies around quasars. 

Here, we present a new technique for studying ab- 
sorbers near luminous quasars, which can be thought of 
as the optically thick analog of the transverse proxim- 
ity effect. Namely, we use background quasar sightlines 
to search for optically thick absorption in the vicinity 
of foreground quasars. Although such projected quasar 
pair sightlines are extremely rare, Hcnnawi ct al. (2006a) 
showed that it is straightforward to select z > 2 projected 
quasar pairs from the imaging and spectroscop y pro- 
vided by the Sloan Digital Sky Survey (SDSS; Yor k et all 
12000(1 . In this work, we combine high signal-to- noise 
rati o (SNR) mode r a.te res olution spectra of the clos- 
est iHennawi et all l)2006ai) projected pairs, obtained 
from Gemini, Keck, and the Multiple Mirror Telescope 
(MMT), with a large sample of wider separation pairs, 
from the SDSS spectroscopic survey, arriving at a to- 
tal of 149 projected pair sightlines in the redshift range 
1.8 < z < 4.0. A systematic search for optically thick ab- 
sorbers in the vicinity of the foreground quasars is con- 
ducted, uncovering 27 new quasar absorber pairs with 
column densities 17.2 < logA^Hi < 20.9 and transverse 
(proper) distances 22 kpc < R < 1.7 Mpc from 
the foreground quasars. 

A handful of quasar-absorber pairs exist in the lit- 
erature, all of which were discovered serendipitously. 
In a study of the statistics of coincidences of opti- 
cally thick absorbers a cross close quasar pair sightlines, 
iD'Odorico et alJ l|2002D discovered one LLS (zabs = 2.12) 
and one DLA (zabs = 2.54) in background quasar spec- 



tra within At; < 1000 km s^ of the foreground quasar 
redshifts, corresponding to transverse proper distances 
of 320 h^^ kpc an d 1.75 M pc, respectively. More 
recently, Adclbcrge r et alJ l|2005(l serendipitously discov- 
ered a faint background quasar (z = 2.92) 49" from a 
luminous (r ~ 16) foreground quasar at z = 2.84, corre- 
sponding to transverse separation R ~ 280 kpc. A 
DLA was detected in the background spectrum at the 
same redshift as the foreground quasar. 

This is the first in a series of four papers on op- 
tically thick absorbers near quasars. In this work, 
we describe the observations and sample selection 
and prese nt 27 new quasar-a bsorber pairs. Paper II 
(Hcnnawi fc Prochaskal l2006al) focuses on the cluster- 
ing of absorbers around foreground quasars and a mea- 
surement of the transverse quasar-absorber correlation 
function is presented. We investigate fluorescent Lya 
emission from our quasar-absorber pairs in Paper III 
l|Hennawi fc Pr ochaska 2006b). Echelle spectra of sev- 
eral of the quasar-LLS systems published here are ana- 
lyzed in Paper IV (Prochaska fc Hennawi 200(3). 

Quasar pair selection and details of the observations 
are described in § El The selection techniques and the 
sample are presented in § |21 A detailed discussion of 
how the systemic redshifts of the foreground quasars were 
estimated is given in § 01 The individual members of the 
sample are discussed in § Cosmological applications 
of quasar-absorber pairs are mentioned in § El a-nd we 
summarize in § |3 

Throughout this paper wc use the best fit WMAP 
(only) cosmological model of Spcrgcl ct al. (2003), with 
flrn = 0.270, flA = 0.73, h = 0.72. Unless otherwise 
specified, all distances are proper. It is helpful to re- 
member that in the chosen cosmology, at a redshift of 
z = 2.5, an angular separation of A9 = 1" corresponds 
to a proper transverse separation of i? = 6 h^^ kpc, and 
a velocity difference of 1500 km s~^ corresponds to a ra- 
dial redshift space distance of s = 4.3 Mpc. For a 
quasar at z = 2.5, with an SDSS magnitude of r = 19, 
the flux of ionizing photons is 130 times higher than 
the ambient extragalactic UV background at an angu- 
lar separation of 60" {R = 340 kpc). Finally, we 
use term optically thick absorbers and LLSs interchange- 
ably, both referring to quasar absorption line systems 
with logTYni > 17.2, making them optically thick at the 
Lyman limit (tll ^ !)■ 

2. QUASAR PAIR OBSERVATIONS 

Finding optically thick absorbers near quasars requires 
spectra of projected pairs of quasars at different redshifts, 
both with z > 2, so that Lya is above the atmospheric 
cutoff. In this section we describe the spectra of pro- 
jected quasar pairs from the SDSS and 2QZ spectroscopic 
surveys as well our subsequent quasar pair observations 
from Keck, Gemini, and the MMT. 

2.1. The SDSS Spectroscopic Quasar Sample 

The Sloan Digital Sky Survey uses a de dicated 2.5m 
telescope an d a large format CCD camera (|Gunn et alJ 
119981 I2006|) at the Apache Point Observatory in New 
Mexico to obtain images in flve broad bands (u, r, i 
and z, cente red at 3551, 4686, 6166, 7480 and 8932 A 
respectively: iFukugita et alJl996l:ISt'oughton et alJ2002() 
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of high Galactic latitude sky in the Northern Galac- 
tic Cap. The ima ging data are pr ocessed by the as- 
trometric pipeline ('Pier et ^aD l2003() and photometric 
pipeline ( LuDton ct al. 2001^7 and a re photometrically 
calibrated to a st andard star network ijSmith et al.ll20n2: 
iHogg et al.ll200lD . Additi onal details on th e SPS S data 
produ cts can be found in lAbazaiian et al.l l)2003l 12004 
|2003). 

Based on this imaging data, spectroscopic targets cho- 
sen by various selection algorithms (i.e. quasars, galax- 
ies, stars, serendipity) are observed with two double spec- 
trographs producing spectra covering 3800-9200 A with 
a spectral resolution ranging from 1800 to 2100 (FWHM 
~ 150 — 170 km s~^). Details of t he spectroscopi c ob- 
servations can b e fouii d in ICastan der et alJ l)2001|) and 
iStoughton et al.l l|2002D . A discussion o f quasar target 
selection is presented in iRichards et alJ l)2002aD . The 
blue cutoff of the SDSS spectrograph imposes a lower 
redshift cutoff of z « 2.2 for detecting the Lya transi- 
tion. The Third Data Release Q uasar Catalog contains 
46,420 quasars jSchneider et al]l2005|) . of which 6,635 
have z > 2.2. We use a larger sample of quasars which 
also includes non-public data: our parent quasar sample 
comprises 11,742 quasars with z > 2.2. Note also that we 
have used the Princeton/MIT spectroscopic reductions"'^'^ 
which differ slightly from the official SDSS data release. 

The SDSS spectroscopic survey selects against close 
pairs of quasars because of fiber collisions. The finite size 
of optical fibers implies only one quasar in a pair with 
separation < 55" can be observed spectroscopically on a 
given platc"'^^. Thus for sub-arcminute separations, addi- 
tional spectroscopy is required both to discover compan- 
ions around quasars and to obtain spectra of sufficient 
quality to search for absorption line systems. For wider 
separations, projected quasar pairs can be found directly 
in the spectroscopic quasar catalog. 

2.2. The 2QZ Quasar Sample 

The 2dF Quasar Redshift Survey (2QZ) is a homo- 
geneous spectroscopic catal og of 44,576 stella r objects 
with 18.25 < bj < 20.85 ijCroom et al.ll200^ . Selec- 
tion of quasar candidates is based on broad band col- 
ors (ubjr) from automated plate measurements of the 
United Kingdom Schmidt Telescope photographic plates. 
Spectroscopic observations were carried out with the 2dF 
instrument, which is a multi-object spectrograph at the 
Anglo- Australian Telescope. The 2QZ covers a total area 
of 721.6 deg^ arranged in two 75° x 5° strips across the 
South Galactic Cap (SGP strip), centered on S = —30°, 
and North Galactic Cap (NGP strip, or equatorial strip), 
centered at 5 = 0°. The NGP overlaps the SDSS foot- 
print, corresponding to roughly half of the 2QZ area. By 
combining the SDSS quasar catalog with 2QZ quasars 
in the NGP we arrive at a combined sample of 12,933 
quasars with z > 2.2, of which 11,742 are from the SDSS 
and 1,191 from the 2QZ. 

The 2QZ spectroscopic survey is also biased against 
close quasar pairs: their fiber collision limit is 30". The 

Available at http:/ /spectro. princeton.edu' 

An exception to this rule exists for a fraction (~ 30%) of 
the area of the SDSS spectroscopic survey covered by overlapping 
plates. Because the same area of sky was observed spectroscopi- 
cally on more than one occasion, there is no fiber collision limita- 
tion. 



fiber collision limits of both the SDSS and 2QZ can 
be partly circumvented by searching for SDSS-2QZ pro- 
jected quasar pairs in the region where the two surveys 
overlap. 

2.3. Keck, Gemini, and MMT Spectroscopic 
Observations 

Another approach to overcome the fiber collision lim- 
its is to use the SDSS five band photometry to search 
for candidate companion quasar s around known, spec- 
troscopically confirmed quasars. IHennawT "eTaLl l|2006ajl 
used the 3.5m telescope at Apache Point Observatory 
(APO) to spectroscopically confirm a large sample of 
photometrically selected close quasar pair candidates. 
This survey discovered both physically associated, binary 
quasars, as well as projected quasar pairs, and produced 
the largest sample of close pairs in existence. 

We have obtained high signal-to-nois e ratio, moder- 
ate resolution spectra of a subset of the iHennawi et alJ 
(2006a) quasar pairs from Keck, Gemini, and the MMT. 
Thus far, 88 quasars with z > 1.8 have been observed, 
which is the operational lower limit for detecting Lya 
set by the atmospheric cutoff. We primarily targeted 
the closest quasar pairs with small separations below the 
fiber-collision limit {A6 < 55"). In some cases other 
nearby quasars or quasar candidates were also observed 
at wider separations from a known close pair. This was 
most often the case with the Keck observations, where a 
multi-slit configuration was used, such that other nearby 
known quasars or quasar candidates could be simultane- 
ously observed on a single mask. Because some of the 
88 quasars we observed are in triples or quadruples, the 
total number of pairs is greater than 44. About half of 
our pairs targeted consisted of projected pairs of quasars 
(Aw > 2500 km s^^) at different redshifts; the rest were 
physically associated binary quasars. 

This spectroscopy program has several science goals: 
to measure small scale transverse Lya forest correlations, 
to constrain the dark energy den sity of the Universe 
with the A lcock-Paczynski test (| Alcock fc~PaczyiisO 
19791 "McDonald fc Miralda- EsTrnd^ 199j 

Hui. Stebbins. & Buries 1999), and to characterize 
the transverse proximity effect. None of these projected 
pairs were specifically targeted based on the presence 
or absence of an LLS. Thus these projected sightlines 
constitute an unbiased sample for searching for optically 
thick absorbers near foreground quasars. 

For the Keck observations, we used the Low Resolution 
Imaging Spectrograph fLRIS: l0ke et al.lll995j) . in multi- 
slit mode with custom designed slitmasks, which allowed 
placement of slits on other known quasars or quasar can- 
didates in the field. LRIS is a double spectrograph with 
two arms giving simultaneous coverage of the near-UV 
and red. We used the D460 dichroic with the 1200 lines 
mm~^ grism blazed at 3400 A on the blue side, resulting 
in wavelength coverage of « 3300 — 4200 A. The disper- 
sion of this grism is 0.50 A per pixel, giving a resolution 
of FWHM- 125 km s"\ On the red side, we used the 
300 lines mm~^ grating blazed at 5000 A, which covered 
the wavelength range 4700 - 10, 000 A, resulting in 2.4 A 
per pixel dispersion or a FWHM- 500 km s"\ All the 
LLSs discovered in the Keck LRIS data were found in 
the blue side spectra, owing to the low redshift (z — 2) 
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of our Keck targets. We used the longer wavelength cov- 
erage on the red side to aid with the identification of new 
quasars and to determine accurate systemic redshifts (see 
§ 0J . The Keck observations took place during two runs 
on UT 2004 November 7-8 and UT 2005 March 8-9. 

The Gemini data were take n with the Gemin i Multi- 
Object Spectrograph (GMOS; iHook et al.ll2004(l on the 
Gemini North facility. We used the B1200_G5301 grating 
which has 1200 lines mm~^ and is blazed at 5300 A. The 
detector was binned in the spectral direction resulting 
in a pixel size of 0.47 A with the 1" slit, correspond- 
ing to a FWHM~ 125 km s"\ The slit was rotated so 
that both quasars in a pair could be observed simultane- 
ously. The wavelength center depended on the redshift 
of the quasar pair being observed. We typically observed 
z ~ 2.3 quasars with the grating center at 4500 A , giv- 
ing coverage from 3750 — 5225 A , and higher redshift 
z ^ Z pairs centered at 4500 A , covering 4000-5250 A . 
The Gemini CCD has two gaps in the spectral direction, 
corresponding to 9 A at our resolution. The wavelength 
center was thus dithered by 15-50A so as to obtain full 
wavelength coverage in the gaps. The Gemini observa- 
tions were conducted over three classical runs during UT 
2004 April 21-23, UV 2004 November 16-18, and UT 2005 
March 13-16. 

At the MMT on UT 2003 December 28-29, we used 
the blue channel spectrograph with the 832 lines mm~^ 
grating at a second order blaze wavelength of 3900 A. 
The CuS04 red blocking filter was used to block con- 
tamination from first order red light. The resolution was 
0.36 A per pixel, or FWHM=75 km s"\ The grating 
tilt again depended on the quasar pair redshift, but a 
typical wavelength center was 3600 A, givmg coverage 
from 3100-4000 A. 

Exposure times ranged from 1800— 7200s, for the Keck, 
Gemini and MMT observations, depending on the mag- 
nitudes of the targets. The SNR in Lya forest region 
varies considerably, but it is always SNR> 5 per pixel 
for the data we consider here. 

One of the projected pairs in our sample, SDSSJ0239- 
0106, was observed with LRIS-B at low resolution be- 
cause because a damped Lya system was detected in the 
background quasars' SDSS spectrum, coincident with the 
foreground quasar redshift, but the SNR of the SDSS 
spectrum was very low. For this observation (spec- 
trum shown in Figure the 300 lines mm~^ grism 
was used with the D680 dichroic. The spectral cover- 
age was from 3500 — 6800 A giving 1.43 A per pixel or a 
FWHM~ 500 km s^^. The date of this observation was 
UT 2005 December 1. 

Finally, we observed the quasar pair SDSSJ1427- 
0121 with t he DEep Imaging M ulti-Object Spectrograph 
(DEIMOS; iFaber et al.ll2003|) on the Keck II telescope 
on UT 2005 May 5. We observed the quasar for two ex- 
posures of 300s with the long-slit mask (0.75" slit), using 
the 600 lines mm~^ grating centered at 7500 A and the 
gg495 blocking filter. The primary reason for this obser- 
vation was to acquire coverage of the Mg II emission line 
of the foreground quasar so as to estimate an accurate 
systemic redshift (see §0}. 



A large number of projected quasar pairs must be 
searched to find optically thick absorbers near quasars. 
For example, the number of absorbers per unit red- 
shift at ~ 2.5 with column densities logiVuT > 19 i s 
dNjdz ~ 0.5 UPeroux et a] . 1 120051 l?rM^a,ra et a1J l2006\ 
In practice, we search within a velocity window Aw = 
1500 km s~^, because of our uncertainty of the fore- 
ground quasar's systemic redshift (see § 0J. The ran- 
dom probability of finding an absorber in the background 
quasar spectrum within Aw of a foreground quasar is 
~ 2%. Thus of order fifty projected quasar pair sightlines 
are required to find a single quasar-absorber pair, with 
absorbers of lower (higher) column density being more 
(less) abundant. As we will see, LLSs are indeed clustered 
around quasars, so that the probability of finding an ab- 
sorber will increase su bstantially for scales smaller t han 
the correlation length ijHennawi fc Prochaskall2'006aj) . 

Our goal is to construct a statistical sample of ab- 
sorbers near quasars which is nearly complete above 
some column density threshold. One approach would 
be to start with complete samples of absorption line sys- 
tems and qua sars, and simply search for quasar-absorber 
pairs. Indeed. IProchaska et al.l l|2005D recently published 
a large complete sample of DLAs found by searching 
the SDSS spectroscopic quasar catalog. However, by re- 
stricting the quasar sample to be complete in a flux or 
volume limited sense, we would substantially reduce the 
number of usable quasars, and furthermore, this would 
not exploit the large number of close projected sightlines 
provided by our Keck/Gemini/MMT observations. In- 
stead, our approach is to search all projected pair sight- 
lines for which an absorber could be detected in the 
background quasar spectrum, at the redshift of the fore- 
ground quasar. This is of course a question about the 
signal-to -noise ratio and resolution of the background 
q uasar spectrum . 

IProchaska et al. ( 2005) demonstrated that the spectral 
resolution (FWHM ~ 150 km s""'^), signal-to-noise ratio 
(SNR), and wavelength coverage of the SDSS DR3 spec- 
tra are well suited to detecting DLAs at z > 2.2. They 
limited their search to objects with SNR> 4 per pixel 
corresponding to roughly r < 19.5, for which they were 
complete at > 95% for column densities logA^Hi > 20.2. 
Here we must be more aggressive because of the limited 
number of projected sightlines and the paucity of known 
LLSs near quasars. Consequently, our sample may not 
yet achieve such a high completeness and it is more sus- 
ceptible to false-positive detections. 

We begin by finding all unique projected quasar pair 
sightlines which have a comoving transverse separation 
oi R < 5 Mpc^^ at the redshift of the foreground 
quasar. The list of potential quasar pair members in- 
cludes all of the quasars in the SDSS-I-2QZ sample, all 
88 of the quasars for which we have Keck/Gemini/MMT 
spectra, and all the quasa rs confirmed f rom APO fol- 
low up spectra published in|Henna wi et al. (2006a). Any 
known quasar can serve as a foreground quasar, provided 
we are confident of its redshift. Only objects which sat- 
isfied the SNR criterion described below could serve as 
background quasars. The 2QZ spectra and the APO 



3. SAMPLE SELECTION 



A comoving distance limit is imposed rather than a proper one 
because the clustering analysis in liHennawi fc Prochaska.2006a.) is 
carried out in comoving units. 
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TABLE 1 

Optically Thick Absorbers Near Quasars 



Name 


Zbg 


Zfg 




R 


^abs 


1 A 1 
\l^V\ 




log iVni 


9uv 


Redshift 


Fg 


Bg 






\ ) 






/"L-m o"!"^ 

^^Kin s ) 


^^Kin S ) 


(cm j 






inst. 


inst. 


SDSSJ0036+0839 


2.69 


2.569 


154.5 


894 


2.5647 


360 


500 


18.95 ± 0.35 


7 


GUI] 


SDSS 


SDSS 


SDSSJOI27+I507I 


2.60 


1.818 


131.0 


794 


1.8188 


30 


300 


18.6 ±0.3 


3 


Mg 11 


LRIS-R 


LRIS-B 




2.38 


1.818 


51.9 


315 


1.8175 


100 


300 


18.9 ±0.3 


13 


Mg II 


LRIS-R 


LRIS-B 


SDSSJ0225-0739 


2.99 


2.440 


214.0 


1251 


2.4476 


690 


500 


19.55 ±0.2 


5 


GUI] 


SDSS 


SDSS 


SDSSJ0239-01062 


3.14 


2.308 


3.7 


22 


2.3025 


540 


1500 


20.45 ± 0.2 


6369 


C IV 


SDSS 


LRIS-B 


SDSSJ0256+0039 


3.55 


3.387 


179.0 


960 


3.387 


20 


1000 


19.25 ± 0.25 


20 


C IV 


SDSS 


SDSS 


SDSSJ0303-0023 


3.23 


2.718 


217.6 


1240 


2.7243 


500 


1000 


18.95 ± 0.2 


8 


G III] 


SDSS 


SDSS 


SDSSJ0338-0005 


3.05 


2.239 


73.5 


436 


2.2290 


960 


1500 


20.9 ±0.2 


13 


C IV-C III] 


SDSS 


SDSS 


SDSSJ0800+3542 


2.07 


1.983 


23.1 


139 


1.9828 


40 


300 


19.0 ±0.15 


488 


Mg II 


LRIS-R 


LRIS-B 


SDSSJ0814+3250 


2.21 


2.182 


10.3 


61 


2.1792 


280 


1500 


18.8 ±0.2 


1473 


Template 


CMOS 


CMOS 


SDSSJ0833+0813 


3.33 


2.516 


103.4 


601 


2.505 


980 


1000 


19.45 ± 0.3 


18 


C III] 


SDSS 


SDSS 


SDSSJ0852+2637 


3.32 


3.203 


170.9 


931 


3.211 


550 


1500 


19.25 ± 0.4 


13 


C IV 


SDSS 


SDSS 


SDSSJ0902+2841 


3.58 


3.325 


183.0 


986 


3.342 


1200 


500 


> 17.2 


34 


GUI] 


SDSS 


SDSS 


SDSSJ1134+3409 


3.14 


2.291 


209.2 


1237 


2.2879 


320 


500 


19.5 ±0.3 


11 


GUI] 


SDSS 


SDSS 


SDSSJ1152+4517 


2.38 


2.312 


113.4 


669 


2.3158 


370 


500 


19.1 ±0.3 


30 


GUI] 


SDSS 


SDSS 


SDSSJ1204+0221 


2.53 


2.436 


13.3 


78 


2.4402 


370 


1500 


19.7 ±0.15 


625 


Template 


CMOS 


CMOS 


SDSSJ1213+1207 


3.48 


3.411 


137.8 


736 


3.4105 


30 


1500 


19.25 ± 0.3 


39 


Template 


SDSS 


SDSS 


SDSSJ1306+6158 


2.17 


2.111 


16.3 


97 


2.1084 


200 


300 


20.3 ±0.15 


420 


Mg II 


LRIS-R 


LRIS-B 


SDSSJ1312+0002 


2.84 


2.671 


148.5 


850 


2.6688 


200 


500 


20.3 ±0.3 


23 


GUI] 


SDSS 


SDSS 


SDSSJ1426+5002 


2.32 


2.239 


235.6 


1397 


2.2247 


1330 


500 


20.0 ±0.15 


19 


GUI] 


SDSS 


SDSS 


SDSSJ1427-0121 


2.35 


2.278 


6.2 


37 


2.2788 


50 


300 


18.85 ±0.25 


7871 


Mg II 


DEIMOS 


CMOS 


SDSSJ1429-0145 


3.40 


2.628 


140.2 


808 


2.6235 


400 


1000 


18.8 ±0.2 


20 


GUI] 


2QZ 


SDSS 


SDSSJ1430-0120 


3.25 


3.102 


200.0 


1100 


3.115 


960 


1500 


20.5 ±0.2 


26 


Template 


SDSS 


SDSS 


SDSSJ1545+5112 


2.45 


2.240 


97.6 


579 


2.243 


320 


500 


19.45 ± 0.3 


30 


GUI] 


SDSS 


SDSS 


SDSSJ1621+3508 


2.04 


1.931 


76.7 


463 


1.9309 


10 


300 


18.7 ±0.2 


12 


Mg II 


LRIS-R 


LRIS-B 


SDSSJ1635+3013 


2.94 


2.493 


91.4 


532 


2.5025 


820 


500 


> 19 


111 


GUI] 


SDSS 


SDSS 


SDSSJ2347+15013 


2.29 


2.157 


47.3 


282 


2.176 


1770 


1000 


> 18.3 


63 


GUI] 


APO 


CMOS 




2.29 


2.171 


223.0 


1329 


2.176 


380 


500 


> 17.2 


8 


Mg II 


SDSS 


CMOS 



Note. — Optically thick absorption line systems near foreground quasars. 
1 In the systems SDSSJ0127+1507 there are two distinct background quasars at ^ = 2.38 and z = 2.60, which show absorption in the vicinity of 
the same foreground quasar at 2 = 1.818. 

^ The foreground quasar for this system has large BAL troughs in the Lya and C IV emission lines. The redshift was computed by comparing 
the peak of G IV, determined by eye, to the shifted wavelength A = 1545.3 A. We apply a conservative redshift uncertainty of ±1500 km . 
^ Voigt profile fits to the Lya absorption in the SDSS spectrum of the background quasar gave log ATjii = 19.55 ± 0.3. An archive echelle spectrum 
of this quasar gives the smaller value which is listed in the table log A^hi = 18.8 ± 0.2. 

^ In the systems SDSSJ2347+1501, there is a single background quasar at z = 2.29 and two foreground quasars at 2 = 2.157 and z = 2.167, 
although the velocity separation is larger than our nominal 1500 km cutoff for the former. 



spc(;tia clo not liave sufBcient resolution or SNR to find 
high column density absorbers, so these quasars could 
only serve as foreground quasars. Furthermore, all of 
our Kcck/Gcmini/MMT spectra easily satisfy our SNR 
criteria, so in practice, we only apply a SNR statistic to 
the SDSS spectra. 

3.1. SNR Statistic 

We define a signal-to-noise statistic SNRbg in the back- 
ground quasar spectrum which is an average of the me- 
dian signal-to-noise ratio blueward and redward of the 

Lya transition at the foreground quasar redshift. 

For the blue side, we begin at the wavelength Abiuc = 
(1 + Xfg)(1215. 67-20) A, and take the median SNR of the 
150 pixels blueward of this wavelength. The 20 A offset 
(4936 km s~^) is applied so that the SNR is not biased 
by the presence of a potential absorber. If there are 
not 150 available pixels blueward of Abiue because of the 
blue cutoff of the spectrum, we take the median of the 
nh\ne > 50 pixels which remain. If less than 50 pixels are 
available, we set SNRtiue = and nuue = 0. 

Similarly, on the red side we begin at Ared = (1 + 



Zfg)(I2I5.67 + 20) A, and take the median SNR of 
the 150 pixels redward of this wavelength. If there 
are not 150 pixels redward of Ared which also have 
A < (1 + Zbg) 1190 A, we compute the median SNR^ed 
of the rired pixels available. Wavelengths larger than 
(l + .2bg)1190 A are avoided because the SNR rises at the 
Lya emission line in the background quasar spectrum. If 
rired < 150, we then also compute the median SNR1275 
of the ni275 = 150 — rircd remaining pixels redward of 
the wavelength A1275 = (1 + Zbg)1275 A, which is free 
of emission lines and a good place to estimate the red 
continuum SNR. Our SNR statistic is defined to be the 
average 

SNRb = '^blue SNRblue -|- rtred SNRygd + ?^1275 SNR1275 
^ ?^blue + '^red + ?^1275 

We require that the foreground quasar's Lya must 
be redward of the background quasar's Lyman limit, 
(1 + Zfg)1215.67 > (1 + Zbg)912, to avoid searching in 
the highly absorbed low SNR region blueward of 912 A. 
A minimum velocity difference between the two quasars 
of Av > 2500 km s~^ is chosen to exclude binary quasars. 
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These pairs with smah velocity separation are excluded 
to avoid confusion about which object is in the back- 
ground and to avoid distinguishing absorption intrinsic 
to the background quasar from absorption associated 
with the foreground quasar. Because the small angu- 
lar separation projected pairs are particularly rare, we 
set a more liberal minimum SNR of SNRbg > 1.5 for 
projected pairs which have (comoving) transverse sep- 
aration R < 1 Mpc. For wider separation pairs 
1 h^^ Mpc < R < 5 Mpc (comoving), we require 
SNRbg > 2. 

3.2. Visual Inspection 

All projected quasar pairs satisfying the aforemen- 
tioned criteria were visually inspected and we searched 
for significant Lya absorption within a velocity window 
of |Au| = 1500 km s~^ about the foreground quasar red- 
shift. This velocity range because it brackets the un- 
certainties of the foreground quasar systemic redshift 
(see Strong broad absorption hne (BAL) quasars 

with large C IV equivalent widths (EWs) were excluded 
from the analyses. Mild BALs were excluded if the BAL 
absorption clearly coincided with the velocity window 
about the foreground quasar redshift which was being 
searched. 

Systems with significant Lya absorption were flagged 
for H I absorption profile fitting. In the SDSS spec- 
tra, all systems which had an absorber with rest equiv- 
alent width W\ > 2 A were fiagged to be fit. We 
adopted a lower threshold of Wx > 1.5 A for the 
Keck/Gemini/MMT spectra, which have higher SNRs 
and slightly better resolution. These equivalent width 
thresholds correspond to column densities of roughly 
logiVni > 19 and logA^ni ^ 18.5, respectively. 

The H I search was complemented by a search 
for metal lines at the foreground quasar redshift, in 
the clean continuum region redward of the Lya for- 
est of the background quasar. The narrow metal 
lines provide a redshift for the absorption line sys- 
tem and, if present, they can help distinguish opti- 
cally thick absorbers from blended Lya forest lines. 
We focused on the strongest l ow-ion transitions com - 
monly observed in DLAs (e.g. iProchaska et alJl2003j) : 
Si II A1260, 1304, 1526, O I A1302, C II A1334, 
Al II A1670, Fe II A1608, 2382, 2600, Mg II A2796, 2803; 
and the strong high-ionization transitions commonly seen 
in LLSs: C IV A1548, 1550 and Si IV A1393, 1402. 
Any systems with secure metal-line absorption were also 
flagged to be fit. 

The Lyman limit at 912 A is redshifted into the SDSS 
spectral coverage for z > 3.2. Although we did not ap- 
ply any specific SNR criteria on the spectra at these 
bluer wavelengths, special attention was paid to pro- 
jected pairs for which the Lyman limit was detectable. 
Systems which showed Lyman limit absorption at the 
redshift of the foreground quasar were also flagged, re- 
gardless of the equivalent width of their Lya absorption 
or the strength or presence of metal lines. 

3.3. Voigt Profile Fitting 

For all of the systems which were flagged by the ini- 
tial visual inspection, we estimated the H I column den- 
sity by fitting the Lya profiles using standard practice. 
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Fig. 1. — Distribution of foreground quasar redshifts, transverse 
separations, and ionizing fluxes probed by the background quasar 
sightlines. The upper plot shows ionizing flux versus proper separa- 
tions, which explains the general trend. The lower plot shows 
foreground quasar redshift versus proper separations and the y-axis 
on the right indicates the wavelength of the Lyo A1215.67 A tran- 
sition at this redshift. The (blue) squares have a Keck (LRIS- 
B) spectrum of the background quasar, (red) triangles have Gem- 
ini (GMOS) background spectra, (magenta) upside down triangles 
have MMT (Blue Channel) background spectra, and (green) circles 
have SDSS background spectra. Filled symbols outlined in black 
have an optically thick absorber at the foreground quasar redshift 
(see Table and open symbols have no absorber. The region to 
the left of the dotted line is excluded by the SDSS fiber collision 
limit of 6 = 55", which explains the paucity of SDSS background 
spectra there. The follow-up Keck/Gemini/MMT spectra probe 
angular separations an order of magnitude smaller than the fiber 
collision limit, allowing us to probe the foreground quasar environ- 
ment down to 20 kpc//i where the ionizing flux is ~ 10,000 times 
the UV background. 



Namely, we over-plotted a Voigt profile on the Lya tran- 
sition, and centered the profile according to the redshift 
of metal-lines, if present. Otherwise, the redshift of the 
absorber was allowed to be a free parameter in the fit. 
The fits are done 'by-eye', which is to say we do not 
minimize a because the error in the fit is dominated 
by systematic uncertainty related to the quasar contin- 
uum placement and line-blending. Conservative error es- 
timates are adopted to account for this uncertainty. In 
all cases, we assume a Doppler par ameter b, which is 
typic al of the high z Lya forest (e.g. iKirkman fc Tvtleil 
2113) ■ In general, the fits are insensitive to the Doppler 
parameter parameter because most of the leverage in the 
fit comes from the damping wings of the line-profile; we 
assume & = 30 km s-^ See Prochaska et al. ( 200^ for 
more discussion on Voigt profile fits to Lya absorption 
profiles. 

The completeness and false positive rate of our sur- 
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vey are sources of concern. Line-blending, in partic- 
ular, can significantly depress the continuum near the 
Lya profile and mimic a damping wing, biasing the col- 
umn density high. To investigate these issues we con- 
sider three objects in our sample for which we have 
both SDSS spectra and independent echelle observa- 
tions. We obtained archived spectra of the background 
quasars SDSSJ 0303-0023, and SDSSJ 1204-^0221, ob- 
served with the High Resolution Echelle Spectrometer 
(HIRES: l\^gt et aUlTool (FWHM w 10 km s"^) on 
the Keck-I telescope, and SDSSJ 1429-0145, observed 
with the Magellan Inamori Kyo cera Echelle spectrograph 
(MIKE: 'Berns tein et alJ 1200 3^ (FWHM « 10 km s"^). 
For SDSSJ0303-0023, the column density of logiVni = 
18.95 ±0.2 measured form the SDSS spectrum is in good 
agreement log A'hi = 18.9 ± 0.2 from the echelle spec- 
trum. Likewise, for SDSSJ 1204+0221 we measured 
logiVm = 19.7 ±0.25 from the SDSS data and logiVm = 
19.7 ± 0.2 from HIRES. However, for SDSSJ1429-0145 
we measure a total column of logA'^Hi — 19.55 ± 0.3 
from the SDSS data; whereas the echelle data gives 
logiVni = 18.8 ± 0.2, a value 2.5cr lower than the SDSS 
value. The source of the error is line blending, but 
we note that this absorber was located blueward of the 
quasars Ly/3 emission line, in a 'crowded' part of the 
spectrum because of the presence of both the Lya and 
Ly/3 forests. 

Based on visually inspecting 149 background quasar 
spectra and and the comparison with the echelle data for 
three systems, we suspect that our survey is ~ 90% com- 
plete for log^Hi > 19.3 for all the Keck/Gemini/MMT 
spectra and the SDSS spectra with SNRbg > 3 (75%), 
or 123 of the 149 spectra we searched. For SDSS spec- 
tra with lower SNRbg > 2, our completeness limit is 
probably closer to logA^Hi > 20. A more careful ex- 
amination of the completeness and false positive rate of 
Lyman limit systems identified in spectra of the reso- 
lution and SNR used here is definitely warranted. Sta- 
tistical studies based on our sample which attempt to 
quantify the abundance of absorbers near foreground 
quasars ( H ennawi fc Prochaska 2006a) , will suffer from a 
'Malmquist' type bias because line-blending biases lower 
column densities upward, and the line density of ab- 
sorbers dN/dz, is a steep function of column density 
limit. Finally, our completeness is likely to be higher and 
the false positive rate lower at 2: ~ 2, as compared higher 
redshifts (z ~ 3), because the mean flux decrement of the 
Lya forest absorption decreases with decreasing redshift, 
and thus line blending causes less confusion. 

3.4. Quasar-Absorber Sample 

We present relevant quantities for the quasar-absorber 
pairs discovered in our survey in Table ^ All systems 
with log A'hi > 17.2 are included; however those systems 
with column densities logA'ni < 18.3 are included as 
LLSs based on the identification of definite Lyman limit 
absorption (i.e. Zfg > 3.2). In the next section (§ 0J, 
we describe in detail how we estimated the foreground 
quasar redshifts and redshift errors which are listed in 
Table ^ The emission line which was used to compute 
this redshift is also noted in the Table. 

The quantity guv = 1 + -Fqso/-Puvb in Tabled is the 
maximum enhancement of the quasars ionizing photon 



flux over that of the extragalactic ionizing background^^, 
at the location of the background quasar sightline, as- 
suming that the quasar emission is isotropic. It is the 
maximal enhancement because we assumed the distance 
is given by the transverse component alone. A discussion 
of how guY was computed is provided in Appendix A. 

The distribution of foreground quasar redshifts, trans- 
verse separations, and ionizing fluxes probed by all of our 
projected pair sightlines is illustrated by the scatter plot 
in Figure^ The filled symbols outlined in black indicate 
the sightlines which have an optically thick absorption 
line system (see Table ^ and open symbols are sight- 
lines without an absorber. Note that the closest pairs are 
predominantly at redshift z ~ 2, becaus e this is where 
quasar pair selection is most efficient fsee lRichards et all 
i2QD2a; Hennawi et al. 2006a). 

A list of tentative absorbers near quasars, for which we 
co uld n ot be sure that logiVni > 17.2, is published in Ta- 
blc lB2l of Appendix B. Higher SNR and higher resolution 
spectra are required to make definitive conclusions about 
these systems; they are a valuable set of targets for future 
research. Coordinates and SDSS fi ve ba nd photometry 
of the quasar pairs in Tables ^ and IB2I are provided in 
Tables lUSl and respectively, of Appendix C. 

To summarize, we searched 149 projected quasar pair 
sightlines with transverse proper distances 22 h^^ kpc < 
R < 1.7 Mpc away from foreground quasars in the 
redshift range 1.8 < z < 4.0. Keck spectra accounted for 
25 of the background quasar spectra, 19 came from Gem- 
ini, 2 from the MMT and the remaining 103 were from 
the SDSS. Of these sightlines, 25 had angular separations 
below the SDSS fiber collision limit (AO < 55") with the 
Keck/Gemini/MMT accounting for all but three of these. 
We discovered 27 LLSs within Av = 1500 km s^^ of the 
foreground quasar redshifts, of which 17 are super-LLSs 
with log A^Hi > 19. 

4. ESTIMATING SYSTEMIC REDSHIFTS 

The primary rest-frame ultraviolet quasar emission 
lines which are redshifted into the optical for z > 2 
quasars are: Lya A1216, N V A1240, C IV A1549, 
C III] A1908, and Mg II A2798. The redshifts determined 
from these lines can differ by up to ~ 3000 km s~^ from 
systemic, due to outflowing/inflowing mater ial in the 
broad line regions of quasars (Gaskell 1982; Tv tler fc FanI 
ll992HVanden Berk e t aL^Ol: Richa rds et al.l2002bfl . A 
redshift determined from the narrow (cr < 200 km s"') 
forbidden emission hues [O II] A3727 or [O III] A5007, 
are the best predictors of systemic redshift; but at z > 2, 
measurements of these lines would require spectra cover- 
ing the near infrared. 

For example, redshifts are most unce rtain when es- 
timated from the C IV emission line: iRichards et alJ 
l)2002bl f) found a median blueshift of 824 km s~ from 
Mg II with a dispersion about the median of 511 km s^^, 
but with a tail extending to blueshifts as large as 
3000 km s^^. Furthermore, Mg II has a median shift 
of 97 km s^^ and a dispersion of 269 km s^^ about 
[O III], which is used to define the systemic frame 
l)Richards et al.l r2002b() . The implied average redshift 

We compare to the UV background computed by F. Haardt 
& P. Madau (2006, in preparation) 
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uncertainty between C IV and [O III] (systemic) is thus 
> 600 km s^ . In addition, the C IV blueshift is luminos- 
ity dependent ?, such that simply adding a median offset 
to a C IV line center can bias redshifts and result in larger 
errors. The quasar redshifts computed by the SDSS spec- 
troscopic pipeline are the result of a maximum likelihood 
procedure which involves f itting of multiple emiss ion lines 
simultaneously (see e.g. iStouehton et aljr2002(l . a pro- 
cedure which will not result in robust systemic redshift 
estimates. 

In light of these issues, we recompute the systemic 
redshifts of all the foreground quasars published in Ta- 
ble n Note that the spectra of foreground quasars used 
to compute these redshifts come from a variety of instru- 
ments (see Table ^ and they have varying SNRs and 
spectral coverage. Some are from the SDSS, for oth- 
ers we have high SNR Keck LRIS-R spectra covering 
A = 5000 - 9500 A, and for others we have only Gemini 
GMOS spectra with 1500 A of coverage centered on 
the Lya forest and typically extending to ~ 1350 A in 
the quasar rest-frame. 

For our spectra with full wavelength coverage (SDSS 
and LRIS-R), we begin by fitting the sum of a Gaussian 
plus a linear continuum (both in log A) to the primary 
emission lines C IV, C III], and Mg II if present, where 
the line-widths and centers are free parameters. We also 
include components for the weaker lines He II A1640 and 
Al IIIA1857, since they can be significantly blended with 
C IV and C III], respectively, thus contributing to the 
background and influencing the placement of the contin- 
uum. Guided by these fits, we calculate the mode of each 
line using the relation mode = 3 x median - 2x mean, 
applied to the upper « 60% of the emission line. This 
is a more robust estimator than the centroid or median 
for slightly skewed profiles in noisy data. Specifically, we 
compute the mode of all spectral pixels within ±1.5cr, of 
the Gaussian line center which have flux 

> '^■^^^ -t- Ca + V p-('°giQ ^-l°gio ^.)V2Tf 

(2) 

where Ai, logj^g-^ij ^^'^ "'ij the amplitude, central 
wavelength, and dispersion, of the best fit Gaussian to 
the ith emission line and C\ is the linear continuum. 
The sum represents the effective background due to other 
nearby lines. 

Given these line centers, a redshift is computed for 
each line. Redshifts are computed using the average ve- 
locity shift from systemic, defined with respect to the 
[O III] eniission lines, where we use t he shifts mea- 
sured b v IVanden Berk et alJ lj2001f) and IRichards et aTl 
Il2002bt) . A median SNR > 5 over the line is required for 
it to be used as a redshift estimator. At high SNR, the 
intrinsic velocity shifts and line asymmetries will domi- 
nate over errors in line centering. In this case, our error 
es timates are mot ivated by the dispersion measurements 
of lRichards et alJ |2002b} . For lower SNR, line centering 
errors can become significant, however we do not explic- 
itly estimate this error contribution. Instead, we simply 
flag low SNR redshift determinations and dilate their er- 
rors by 'hand'. Redshifts and errors are assigned to the 
27 foreground quasars in Tabled according to the follow- 
ing procedure: 



- Use Mg II if it is present {z < 2.2). The redshift 
error is ±300 km s^^ (6 systems). 

- If Mg II is not present or has low SNR, C III] is 
used. The error is ±500 km s^^^^ (13 systems) 

- If neither Mg II nor C III] can be used, the redshift 
is computed from C IVThe error is ±1000 km s^^^^ 
(3 systems). 

- If neither Mg II, C III], or C IV have SNR < 5, aU 
are fit simultaneously for the redshift, using Gaus- 
sians plus continua (similar to SDSS pipeline pro- 
cedure). The error is ±1500 km s"^ (1 system). 

- If C IV is unavailable because of limited spectral 
coverage (i.e. Gemini GMOS spectra), the redshift 
is computed from cros s correlation with a com - 
posite quasar template ijVanden Berk et ani2001(l . 
The redshift error is ±1500 km s~^ (4 systems). 

5. NOTES ON INDIVIDUAL ABSORPTION SYSTEMS 

The Keck, Gemini, and SDSS spectra of both the fore- 
ground and background quasars from Tableware shown 
in Figure 121 The right panels show show closeups of the 
H I absorption in the background quasar spectra and our 
best fit Lya profiles. Below we describe how each system 
was fit and highlight interesting details. 
SDSSJ 0036+0839 The SDSS spectrum of the back- 
ground quasar reveals strong absorption features of 
C II A1334, Si IV A1393,1403, C IV A1548, 1550, 
Si II A1526 and possibly O I A1302 and Si II A1304. 
A strong feature is coincident with Si II A1260, but this 
could be a Lya forest absorption line. Centered on the 
redshift of the metal lines at ^ = 2.5647 ± 0.00005, 
our Voigt profile fit to the H I Lya gives logiVni = 
18.95 ± 0.35. A redshift of Zfg — 2.569 is measured 
from the foreground quasars C III] emission line with 
an error of ±500 km s^^. The velocity difference of 
]Aw] — 360 km s~^ from the absorber is within this error 
estimate. 

SDSSJ 0127+1507 While most of the systems in Ta- 
ble n contain only two quasars, there are five quasars 
within 5' of each other in the region surrounding 
SDSSJ 0127+1507 at redshifts 2.60, 2.38, 2.18,2.08, 
and 1.81 (denoted ABCDE). Along these sightlines, we 
identify two absorption systems in SDSSJ 0127+1507A 
(zbg = 2.60; A0 = 131.0") and SDSSJ 0127+1507B 
(zbg ~ 2.38; A9 = 51.9") with Lya profiles indicating 
an optically thick absorber near the redshift of the fore- 
ground quasar SDSSJ 0127+1507E (zfg = 1.818). No sig- 
nificant metal-line absorption is observed in either back- 
ground quasar spectrum. Voigt profile fits to the two 
profiles give z = 1.819 ± 0.001 and logTVm = 18.6 ± 0.4 
toward SDSSJ0127+1508A and z = 1.1818 ± 0.001 with 
logTVm = 18.9 ±0.3 toward SDSSJ0127+1507B. We fur- 
ther identify additional strong Lya absorption features 
in the spectrum of the background quasar J0127+1507C 

The dispersion of C III] about systemic has not been quanti- 
fied, but it is hkely comparable to that of Mg II (G. Richards, pri- 
vate communication 2005), motivating our choice of ±500 km s~^. 

A conservative uncertainty of ±1000 km s~^ is assumed for 
C IV to account for redshift errors due to asymmetries, self- 
absorption, BAL features, and large intrinsic velocity shifts. 
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(zbg ~ 2.18; A^? = 35"), but we cannot precisely de- 
termine their iVni values, and hence whet her an op- 
tically thick absorber is present (see Table IB2|I . The 
fourth background quasar SDSSJ 0127-I-1507D ((zbg = 
2.18; M) ~ 163.1") shows no signs of significant Lya ab- 
sorption at the redshift Zfg = 1.818. Coordinates and 
photometry for all five members of this projected group 
of quasars are given in Table lUSl of Appendix B. 
SDSSJ 0225-0739 We detect strong Si II A1526 
and Al II A1670 in the background quasar spectrum, 
which identify the redshift of the absorber at z = 
2.4476 ± 0.0007, coincident with the foreground quasar 
at Zfg = 2.440. Our fit to the Lya profile is compli- 
cated by O VI emission from the quasar but yields the 
value logA'^Hi = 19.55 ± 0.2 dex. The velocity differ- 
ence of I Aw I = 690 km s^^ from the absorber is in ex- 
cess of our estimated foreground quasar redshift error of 
±500 km s"\ 

SDSSJ 0239-0106 The SDSS spectrum of the back- 
ground quasar reveals strong absorption features of 
Si II A1526, Fe II A1608 and Al II A1670 at a redshift 
z = 2.3025 ± 0.0007. Our estimate of the foreground 
quasar redshift is Zfg — 2.308 is particularly uncertain 

(±1500 km s~^) because of large BAL features. There is 
an apparently strong Lya profile at this redshift in the 
SDSS spectrum but the SNRfg = 1.8 in this region is 
too low for a meaningful fit. As described in §12 we ac- 
quired a low-dispersion LRIS-B spectrum of this pair to 
improve the fit, and it is these spectra which are shown 
in Figure 121 The metal-lines are detected in the LRIS-B 
spectrum of the background quasar, and a fit to the Lya 
profile gives log A'hi — 20.45 ± 0.2 dex. This absorber, 
therefore, satisfies the damped Lya threshold (log A^hi > 
20.3). It was not included in the Prochaska ct al. (2005.) 
compilation, however, because of the low SNR of the 
SDSS spectrum in this region. 

SDSSJ 0256+0039 The SDSS spectrum of the back- 
ground quasar shows a weak C IV absorption system at 
z — 3.387±0.001 and a Lyman limit system with consis- 
tent redshift. There may also be weak absorption from 
Si II A1260 and Si IV A1393, 1402. An examination of 
the Lya profile at this redshift shows a significant fea- 
ture, centered at zj^ya — 3.387. A fit indicates an H I 
column density of log A'hi = 19.25 ± 0.25 is allowed by 
the data, but this value is very uncertain because of the 
low SNR and lack of evidence for significant damping 
wings. The presence of an LLS sets a lower limit of to 
the H I column density logA^Hi > 17.2, unless the limit 
corresponds to another absorber nearby in redshift. The 
redshift of the foreground quasar Zfg — 3.387 is consistent 
with the redshift of the absorber. 

SDSSJ 0303-0023 There is an absorption system 
showing very strong Si IV and C IV lines at z = 2.7243 ± 
0.0007 in the spectrum of the background quasar, which 
is 500 km s""'^ away from the foreground quasar at z = 
2.718 (±1000 km s~^). Weak absorption consistent with 
Al II A1670 is observed at this redshift, but no signifi- 
cant detection of Fe II A1608 (Wa < 100mA). A strong 
Lya profile is consistent with the C IV lines, although 
with nominal line centroid blueward of the metal-line 
absorption {zi^ya = 2.723). If we constrain the fit by 
the C IV redshift, the red wing of Lya provides the only 
constraint on the fit and we find logiVni = 18.95 ± 0.20. 



If we allow the redshift to be a free parameter, the best 
fit solution is « 0.1 dex larger. As mentioned in the 
previous section, we have independent echelle data for 
this background quasar and a fit to the Lya profile gave 
log A^Hi — 18.9 ± 0.20, consistent with our measurement 
from the lower resolution SDSS spectrum. 
SDSSJ0338-0005 This damped Lya system shows a se- 
ries of metal- line transitions including the Mg II doublet. 
The redshift of the gas is z = 2.2290 ± 0.0007 and the 
foreground quasar is at z = 2.239, but with large error, 
±1500 km s~^. The a bsorber is a memb er of the SDSS- 
DR3 DLA catalog and lProchaska et all (HflOS) report an 
H I column density of logTVni = 20.9 ± 0.15 dex. 
SDSSJ0800+3542 This absorber is notable for show- 
ing strong N V A1238, 1242 absorption, which is rarely 
ever observed in LLSs. There is also significant absorp- 
tion due to Si II A1260 and C II A1334, as well as the 
Fe II A1608 and Al II A1670 transitions. The metal line 
profiles are centered at z = 1.9828 ± 0.0007 and are ac- 
companied by a Lya profile with logiVni — 19.0±0.15, 
according to our fit. The absorber redshift is in good 
agreement with that of the foreground quasar, Zfg — 
1.983. 

SDSSJ0814+3250 Similar to SDSSJ0800+3524, this 
absorption system shows remarkably strong N V absorp- 
tion in the spectrum of the background quasar, as well 
as several other metal transitions (Si II, C II, Si IV). 
The redshift of these metal-hnes is z = 2.1792 ± 0.0007. 
At the corresponding position of H I Lya there is sig- 
nificant absorption, although oddly, the largest optical 
depth occurs 1.2 A away from the line centroid predicted 
by the metal lines. We have fit the Lya profile with 
two H I components, one centered at z = 2.1792 and 
an additional component at z = 2.1778. Our best so- 
lution has logA'ni — 18.50 ± 0.2 for each component. 
There is significant degeneracy between the two compo- 
nents, especially if one does not demand a component at 
z — 2.1792. However, the total column density is well 
constrained and has a value of logA^Hi = 18.80 ±0.2 dex. 
The foreground quasar is at ztg = 2.182 ± 1500 km s^^, 
which is I Aw I = 280 km s~^ away from the absorber (at 
z = 2.1792), but well within the uncertainty of the fore- 
ground quasar redshift. 

SDSSJ0833+0813 This system shows likely Si II A1526 
and possible C IV absorption a,t z — 2.505 ± 0.001 and a 
strong corresponding Lya profile. If we adopt the metal 
line redshift, our best fit is a multi-component model 
with log A^Hi = 19.45 ± 0.3. But, there is a small chance 
that this profile is the blend of several logA^Hi ~ 18 ab- 
sorbers and that the total logA'ni < 19. The metal line 
redshifts are 980 km s~^ from the foreground quasar at 
Zfg = 2.516, but this velocity difference is comparable to 
its redshift error ±1000 km s~^. 

SDSSJ0852+2637 The foreground quasar redshift 
(z = 3.203) is sufficiently high that the SDSS spectrum 
covers the Lyman limit region. We do detect a break 
due to the Lyman limit at 912 A and its correspond- 
ing Lya profile at z = 3.211±0.001, in the background 
quasar spectrum. There are no detected metal-line tran- 
sitions, which is not uncommon for Lyman limit systems 
at the resolution of the SDSS spectra. We estimate the 
column density to be log A^hi — 19.25 ± 0.4. There is a 
|Au| = 550 km s~^ between the absorber and the fore- 
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ground quasar redshift (zfg = 3.203), but this is within 
our estimated ±1500 km s^^ margin of error on the lat- 
ter. 

SDSSJ0902+2841 Similar to SDSSJ0852+2637, the 
background quasar spectrum covers the Lyman limit at 
the redshift of the foreground quasar. The SDSS spec- 
trum of the background quasar reveals a Lyman limit 
system at 2; = 3.342±0.005, corresponding to an offset by 
|Au| = 1200 km s~^ from the foreground quasar redshift 
Zfg = 3.325. This offset is more than twice the estimated 
redshift error (±500 km s~^) from the C III] emission line 
of the foreground quasar. No significant metal-line ab- 
sorption is identified. There is a Lya absorption profile 
with centroid consistent with the Lyman limit redshift. 
There appears to be residual flux in the line profile, but 
the data has insufficient SNR in this region to say defini- 
tively. While a column density log A'hi > 19 is permitted 
by the data, we suspect the actual value is less than this. 
Thus we set a conservative lower limit of log A^hi > 17.2 
based on the presence of a limit. 

SDSSJ1134+3409 This quasar exhibits a very strong 
C IV system at 2 = 2.2879±0.GG5 and possible detections 
of Si II A1526 and Al II A1670 in the background quasar 
spectrum. There is a corresponding Lya profile with 
rest equivalent width Wx — 4.9 ± 0.2A. Our fit to this 
profile gives log A^hi — 19.5 ± 0.3 if we tie the redshift to 
the C IV absorption. We derive a +0.3 dex larger value 
if we allow the redshift to vary as a free parameter to 
z — 2.2865. Either redshift is in good agreement {\Av\ < 
400 km s^ I) with the foreground quasar at Zfg = 2.291. 
SDSSJ1152+4517 The background quasar spectrum 
has a strong Lya absorption feature aX z — 2.3158 ± 
0.0010 and a C IV absorption system a,t z — 2.31323 ± 
0.0007, significantly offset from the Lya feature. Our 
profile fit to Lya gives log A^hi = 19.1 ± 0.3 assuming the 
redshift z = 2.3158. The absorber is offset by \Av\ = 
370 km s^ from the foreground quasar at Zfg = 2.312 ± 
500 km s^^, but this is within the redshift error. 
SDSSJ1204+0221 For this quasar, we have both 
CMOS and SDSS spectra. Figure □ shows the CMOS 
spectra (left) of both pairs, but we used the SDSS data 
to fit the Lya profile (right) , because this spectrum gave 
coverage of more of the metal lines. A series of strong 
metal transitions were identified at 2 = 2.4402 ± 0.0003 
and our fit to the Lya profile gives log A^hi = 19.7 ±0.25. 
A fit to the column density using an echelle spectrum 
(HIRES) of this object gives the same value. The ab- 
sorber redshift is offset by |Ai;| = 370 km s^^ from that 
of the foreground quasar Zfg — 2.436 ± 1500 km s^^, but 
the error on this determination is large. 
SDSSJ1213+1207 A strong Lya absorption feature 
is identified at z — 3.4110 ± 0.0007 in the background 
quasar spectrum, which is in good agreement with the 
foreground quasar redshift Zfg = 3.411. A likely coin- 
cident C IV absorption feature is also identified. The 
SDSS data provide coverage of the Lyman limit at this 
redshift, and the break is detected, although the SNR of 
the data in this region is poor. Our Voigt profile fit to 
the Lya feature gives log A'hi — 19.25 ± 0.3. 
SDSSJ1306+6158 The LRIS-B spectrum of the back- 
ground quasar shows a damped Lya system at z = 
2.1084 ± 0.0007, consistent (|Ai;| = 200 km s~^) with 



the redshift of the foreground quasar, Zfg — 2.111 ± 
300 km s~\ A series of metal-line transitions outside 
the Lya forest are identified, which constrain the red- 
shift well. Our single component fit to the Lya profile 
gives log A^Hi = 20.30 ± 0.15. 

SDSSJ1312+0002 The low SNR spectrum of the 
background quasar shows the following metal-line tran- 
sitions at z = 2.6688 ± 0.0007: O I A1302, 
Si II A1304, C II A1334, Si IV A1393, 1402, Si II A1526, 
C IV A1548, 1550, Al II A1670, and likely Fe II A2382. 
A strong Lya absorption profile is consistent with this 
redshift. Although the SNR is low (sa 2 per pixel), the 
H I column density is consistent with the DLA threshold 
log A^Hi = 20.3 ± 0.3. This pro file was fit with a similar 
value bv lProchaska et aP l|2005t ) but not included in their 
sample because of the low SNR. The absorber redshift is 
offset by |Ai;| — 200 km s~^ from that of the foreground 
quasar Zfg = 2.671, but within its error (±500 km s"'). 
SDSSJ1426+5002 We identify an absorption system in 
this SDSS background quasar spectrum with very strong 
metal- lines including detections of the weak Si II A 1808 
and Zn II A2026 tran sitions at z = 2.2247 ± 0.0007. 
'Proch aska et alJ l)2005(l reported an A'hi column density 
for this absorber below the DLA threshold (logA^ni > 
20.3) and we similarly find log A^hi = 20.00 ± 0.15 here. 
This low column density is rather remarkable relative to 
the the strong metal-line absorption. The redshift of the 
absorber is offset by |Ati| — 1330 km s^^ from the fore- 
ground quasar redshift ztg — 2.239, which is more than 
twice our estimated redshift error of ±500 from the C III] 
emission line of the foreground quasar. 
SDSSJ1427-0121 Our high SNR CMOS spectrum of 
the background quasar reveals several metal-line absorp- 
tion features redward of Lya. Near the redshift of the 
foreground quasar, we identify Si II A1260, C II A1334, 
and Si IV absorption. These lines show multiple com- 
ponents spanning over 1000 km s~^. We identify similar 
structure in the Lya profile, although, the bluest compo- 
nent of the metal lines does not correspond to significant 
H I absorption, similar to the case for SDSSJ 0814+3250. 
We have fit the Lya complex with 3 components at z — 
2.2762, 2.2788, and 2.2828 with log AThi = 14, 18.85, and 
18.80, respectively. The low column density log A'hi = 14 
component results from the lack of significant optical 
depth in the bluest component. Because the features 
are well separated, there is minimal degeneracy in the 
fit, although the uncertainty is still « 0.25 dex for each 
component. The two features at z = 2.2788 and z = 
2.2828 differ in redshift by 366 km s~^, greater than the 
300 km s^^ maximum value suggested bv lO'Meara et all 
(|2006') in combining multiple component absorption sys- 
tems to one. Thus, we take the component at z = 2.2788 
with logA'ni — 18.85 as the quasar companion, since 
it is closest to the foreground quasar systemic redshift 
Zfg = 2.278, differing by only \Av\ = 50 km s"\ 
SDSSJ1429-0145 The SDSS spectrum of the back- 
ground quasar shows an absorption system with large 
EW metal-line features of Si II A1526, C IV A1548, 1550, 
Fe II A1608 and Al II A1670 at a redshift z = 2.6235 ± 
0.0007, in relatively good agreement with the foreground 
quasar redshift at Zfg = 2.628. We identify a strong Lya 
profile with rest equivalent width W\ = 3.8A; however, 
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the centroid of this profile hes blueward of the peak in the 
optical depth of the metal-line transitions. There is also 
evidence for weak absorption in the metal-line profiles 
redward of the dominant absorption. Fitting the profile 
with two components, we derive logTYni = 19.3 ± 0.3 at 
z = 2.6235 and log A^m = 19.2 ± 0.3 at z = 2.6267. Sim- 
ilar to the absorbers associated with SDSSJ0814-I-3250, 
the two values for the two components are degenerate, 
but their total H I column density, A^hi = 19.55 ± 0.3 
is well constrained. As mentioned in § |31 we have 
also examined echelle spectra of this background quasar, 
and a significantly lower column density was measured, 
logA^Hi = 18.8 ±0.2, with the difference being caused by 
unresolved line-blending in the SDSS spectrum. It is this 
lower value which is listed in Tabled 
SDSSJ1430-0120 The background q uasar exhibits a 
pair o f damped Lya systems from the iProchaska et al.l 
ij2005H compilation at z = 3.0238 and z = 3.115 with 
log AThi = 21.10 ± 0.2 and log A^m = 20.50 ± 0.2, respec- 
tively. The DLA at z = 3.115 is |Au| = 960 km s"^ 
away from the foreground quasar Zi^ = 3.102; however 
the foreground quasar redshift has a large associated er- 
ror ±1500 km s~^. Also, there is no apparent metal- 
line absorption for the DLA coincident with the quasar, 
and the redshift is not tightly constrained by the Lya 
profile fit: we estimate an uncertainty of bz — ±0.002 
(±150 km s"^). 

SDSSJ1545+5112 We identify a strong Lya profile at 
z — 2.243 ± 0.001 in the background quasar spectrum of 
this pair, in good agreement with the foreground quasar 
redshift Zfg = 2.240. Weak absorption features consistent 
with O I A1302, C II A1334 and possibly C IV are present; 
but we note the absence of any obvious Si II A1260, call- 
ing these identifications into question, because the oscil- 
lator strengths of these transitions are comparable. A 
Voigt profile fit to the Lya gives log A^hi = 19.45 ± 0.3. 
SDSSJ1621+3508 Our LRIS-R spectrum of this bright 
background quasar exhibits a Mg II system at z = 
1.9309 ± 0.0007, in good agreement with the foreground 
quasar redshift Zfg = 1.931. It also shows metal ab- 
sorption from C II, C IV, and Si II. A correspond- 
ing strong H I Lya absorption profile is identified in 
our LRIS-B spectra and the best fit column density is 
logA^Hi = 18.7 ±0.2. 

SDSSJ1635+3013 The SDSS spectrum of this 
background quasar shows metal-line absorption from 
Si II A1526 and Al II A1670 at a z = 2.5025 ± 0.001. 
The foreground quasar at Z{g = 2.493 ± 500 km s~^ is 
offset from the absorber by |A?;| = 820 km s~^ which 
is significant compared to the foreground quasar redshift 
error. The Lya profile at the redshift of the metals lines 
appears significantly blended with other coincident Lya 
absorption transitions. The H I value is thus poorly con- 
strained. We obtain a fit of logA^ni = 19.35 ± 0.3 but 
report a conservative lower limit of log Ahi > 19. 
SDSSJ2347+1501 This triple system consists of a 
background quasar at SDSSJ2347+1501A at Zbg = 2.29, 
and two foreground quasars: SDSSJ2347+1501B (zfg 
2.157; AO = 47.3) and SDSSJ2347+1501C (zfg = 2.171; 
Ae = 223.0). We identify S IV absorption at z = 2.176± 
0.001, near the redshifts of the foreground quasars. There 
is a corresponding strong Lya absorption feature, but 
it appears to be a blend of two components. Our best 



fit profile, which is not very well constrained, gives 
logA^Hi = 18.55 at z = 2.174 and logA^m = 18.1 at z = 
2.177. These values are highly correlated but the absence 
of significant damping wings limits the total H I column 
density to be log A^hi < 19, and we place a lower limit on 
the total column density of log A'hi < 18.3. The absorber 
is offset by |At;| = 1770 km s"^ from SDSSJ2347+1501B 
at Zfg = 2.157 ± 1000 km s"\ and |Aw| = 380 km s"^ 
from SDSSJ2347+1501C at Zfg = 2.171 ± 500 km s"\ 
The large offset from SDSSJ2347+1501B exceeds the 
limit of our search (±1500 km s^^). If part of this off- 
set is due to Hubble flow, the quasar would be separated 
from the absorber by >> 5h~^ Mpc in the radial direc- 
tion, much larger than the R = 282 kpc transverse 
separation. 

6. COSMOLOGICAL APPLICATIONS OF OPTICALLY 
THICK ABSORBERS NEAR QUASARS 

The cosmological applications of optically thick ab- 
sorption line systems near quasars are many. They in- 
clude measurement of the clustering of absorbers around 
quasars, determining the spatial distribution of neutral 
gas in LLSs and DLAs, fluorescent Lya emission from 
LLSs, constraints on the emission geometry and lifetimes 
of quasars, and studying the enrichment and velocity 
fields in the environs of luminous quasars. We discuss 
each of these in turn. 

6.1. Clustering of Absorbers around Quasars 

A conspicuous feature of Figure ^ is the high cover- 
ing factor of absorbers for small separations. Six out of 
eight projected sightlines with R < 150 h^^ kpc have a 
LLS coincident with the foreground quasar, of which four 
are super-LLSs (log A^hi > 19). We previously estimated 
that the probability of a random quasar-super-LLS coin- 
cidence is ~ 2% (see beginning of §|31l; thus the average 
number expected in eight sightlines is ~ 0.16. This fac- 
tor of ~ 25 increase in the number of small scale quasar- 
absorber coincidences provides significant evidence that 
these absorbers are strongly clustered around quasars. 
It is particularly interesting that in the projected pairs 
which show absorption, an LLS is seen in the transverse 
direction only - none of the foreground quasars show an 
intrinsic or proximate optically thick absorber along the 
line of sight. 

Although they are rare, proximate absorbers with 
Zahs ~ are occasio nally observed in single lines of 
sight. In a recent studv. iRiissell et iiD l)20051 found that 
the abundance of DLAs is enhanced by a factor of ~ 2 
near quasars, which they attributed to the clustering of 
the DLA counterpart galaxies around quasars. The high 
covering factor for R < 150 kpc in Figure sug- 
gests a significantly larger enhancement in the transverse 
direction. We measure the transverse quasar-absorber 
correlation function and compare it to the abundance of 
proximate absorbers in the second paper of this series 
ijHennawi fc Prochask£3l2006a|) . 

An asymmetry between the strength of cluster- 
ing in the radial and transverse direction could be 
the hallmark of anisotropy or obscuration of quasar 
emission. If the emission were highly anisotropic, 
optically thick absorbers along the linc-of-sight might 
be evaporated by the ionizing flux; whereas, trans- 
verse absorbers could lie in shadowed regions, and 
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hence survive. This speculation gains some credibihty 
in Hght of the recent null detections of the trans- 
verse proximity effect i n the L y g f orests of projected 
quasa r p airs llCrottsI 119891: IPo brzvcki & Bcchtold 
iFernand cz-Soto. Barcon s. Carballo. fc Webhl 

T99a .Liskc fc WiUigen mml 

Schirber. Miralda-Escude. fc McDonald! I200I ICrofT 



2001 " Dut see Jakobsen et al.2003). Although these stud- 
ies are each based only on a handful of projected pairs, 
they all come to similar conclusions: the amount of (op- 
tically thin) Lya forest absorption, in the background 
quasar sightline near the redshift of the foreground 
quasar, is larger than average rather than smaller - 
the opposite of what is expected from the transverse 
proximity effect. Our Keck, Gemini, MMT, and SDSS 
spectra will be used to stud y the transverse prox imity 
effect in an upcoming paper l)Hennawi et alJl2b061j) . 

Finally, if obscuration of the quasars emission is in- 
deed the explanation for the anisotropic clustering pat- 
tern suggested by Figure ^ we would naively expect the 
transverse covering factor to be approximately equal to 
the average fraction of the solid angle obscured. Stud- 
ies of Type II quasars and the X-ray background sug- 
gest that quasars with luminosities comparable to our 
foreground quasar sample (Mq < —23) have ~ 30% of 
the s o lid angle obscured (Ueda et al. 2003; Barec r'e't al.l 
120051 iTreister fc Urrvl 120051) . althouerh these estimates 
are highly uncertain. This seems to be at odds with the 
high covering factor (6/8) for having an absorber with 
logiVni > 17.2) that we observe on the smallest scales 
{R < 150 h^^ kpc), although our statistics are clearly 
very poor. 

6.2. Fluorescent Lya Emission: Shedding Light On 
Lyman Limit Systems 

Are there LLSs or DLAs that can self-shield the in- 
tense ionizing flux of a nearby quasar? How do we know 
if these systems are being illuminated? - By observ- 
ing the fluorescent glow of their recombination radia- 
tion. Optically thick H I clouds in ionization equilib- 
rium with a photoionizing background, re-emit ~ 60% 
of the ionizing radiation they absorb as fluorescent Lyg 



recombination line photons llGould fc Weinberg 



Zheng fc Mira,1da,-Fscud_el I2002t ICa,nta,limo et al 



Adelberger et alJl2005t iKollmeier et all 120061) . All at- 
tempts to detect this fluorescent radiation from LLSs 
illuminated b y the ambient UV background have been 
unsuc cessful ("Bun ker et all 119981: iMarleau etHl Il999t 
iBecke r 2005), despite > 60 hour integrations on 10m tele- 
scopes (^Beckeiii2005j) . However, if a nearby quasar illu- 
minates an absorber, as could be the case for the quasar- 
absorber pairs in Tabled the fluorescent surface bright- 
ness is enhanced by a factor of girv, making it much eas- 
ier to detect. Indeed, lAdelberger et al.l l)2005j) reported 
a possible detection of fluorescence from a quasar-DLA 
pair at z = 2.84, which was serendip itously discovered 
in the ir LBG survev l|Steidel et al.ll999H Adelberger et al.l 
EqqI. The ionizing flux is enhanced by a factor of guv — 
4850, owing to their extre mely luminous foreground 
quasa r r ~ 16 (but also see iFrancis fc Bland-HawthornI 
120041 who failed to detect fluorescence around a bright 
B = 17.6 quasar). 

In Table^we publish seven new quasar-absorber pairs 
with enhancements guy > 100, with largest being ^ 7900 



for SDSSJ1427-0121. This corresponds to fluorescent 
surface brightnesses of fj^hya = 19.5 — 24.3 mag arcsec"^, 
if the foreground quasars emit isotropically. This po- 
tential 5-10 magnitude increase in the expected sur- 
face brightness would constitute a breakthrough for the 
study of fluorescence from optically thick absorbers. It 
is intriguing that two of the three systems which have 
enhancements of guv > 1000 (SDSSJ0814-h3250 and 
SDSSJ1427-0121) for which we have high SNR ratio 
moderate resolution spectra, both have H I Lya ab- 
sorption misaligned with the regions of highest metal 
line optical depth. These peculiar absorption profiles 
are suggestive of emission in the LLS trough (see Fig- 
ure |21 and §|31). This suggestion is particularly intriguing 
when one considers that proximate DLAs observed along 
the line of sight to single quasars, seem to preferentially 
exhibit Lya emission superimposed on the Lya ab sorp- 
tion trough iMoller et allligM lEllison et al.ll200l . We 
present a detailed discussion of fluorescent emission from 
our quasar-LLS pairs in the third paper of this series 
enn awj fc Pr ochaska 2006b). 

Finally, because the fluorescent emission from a LLS 
comes from the a thin i onized 'skin' where tll ^ 1 
lIGould fc Weinberg! IT99I . the size of this region can 
be measured if fluorescent emission is detected. Indeed, 
Edelberger et al. (2005) constrained the emitting region 
from their fluorescing DLA to be > 1.5", or a proper size 
of ^DLA = 8.5 kpc. A statistical study of fluorescent 
emission from LLSs illuminated by quasars, can measure 
the distribution of sizes for optically thick absorbers, as 
well as the relationship between size and observed H I 
column. The fluorescent flux also measures the ionizing 
flux at the self-shielding boundary, which, when com- 
bined with the size, constrains the density and distribu- 
tion of gas in LLSs and DLAs, providing physical insights 
into the morphologies of high redshift galaxies which are 
difficult to derive from the resonance line observations of 
absorption line spectroscopy. 

6.3. Constraining Quasar Lifetimes 

In addition to serving as a valuable laboratory for 
studying fluorescent emission and the size and nature of 
LLSs, the detection of fluorescence from an LLS illumi- 
nated by a quasar can place interesting constraints on the 
lifetimes of individual quasars. A LLS near a quasar acts 
like a mirror, 'reflecting' ~ 60% of the ionizing photons 
absorbed from the quasar towards us. However, this re- 
flection can only be detected provided that ionizing flux 
of the quasar has not changed significantly for a time 
longer than the light crossing time to the LLS. The light 
crossing time corresponding to 60" (1.2 h~'^ Mpc) from 
a z = 2.5 quasar is 1.1 x 10^ yr. 

Currently, the lower limit on the intermittency of 
quasar em ission comes from observations of the proxim- 
ity effect l|Baitlik et al.lll988l IScott et al.ll200nD . in the 
Lya forests near quasars. The presence of a proxim- 
ity effect implies that the IGM has had time to reach 
ionization equilibrium with the quasars increas ed ioniz- 
ing fl ux, giving the lower limit iburst ^ 10* yr ({Martini 
|2004|) . However, owing to the large amplitude of intrinsic 
Lya forest fluctuations, many quasars must be averaged 
over to detect the proximity effect, hence the intermit- 
tency of indivi dual quasars cann ot be constrained with 
this technique (jScott et alJl2000tl . A positive detection 



QUASARS PROBING QUASARS 



13 



of fluorescence allows one to place much stronger con- 
straints on quasar intermittency and lifetimes, with the 
added advantage that the lifetimes of individual quasars 
can be constrained. 

6.4. Probing the AGN Environment with 
High- Resolution Spectroscopy 

If the background quasar is sufficiently bright {R < 
19), then one can obtain a high-resolution spectrum on 
a 10m class telescope. These data provide more accurate 
measurements of A'hi for absorbers with logA^ni < 20. 
More importantly, these data yield accurate ionic column 
density measurements critical to the analysis of physical 
conditions (e.g. ionization state, gas density, metallicity). 
It is noteworthy that several of the absorbers exhibit 
the N V doublet in the low-dispersion data. This ion 
is rarely observed in quasar absorption line systems and 
suggests an environment with a significantly enhanced 
UV radiation field. The spectra also reveal the velocity 
field of gas in the vicinity of bright AGN where feedback 
processes (e.g. jets, merger activity) may be important. 
In the case of SDSS1204 we measure a velocity width 
Av « 600 km s"^ for the Si+, C+, and ions. We 
will present and analyze these observatio ns in Paper IV 
of the series ijProchaska & Hennawill2'00(iD . 

7. SUMMARY 

We searched 149 moderate resolution background 
quasar spectra, from Gemini, Keck, MMT, and the 
SDSS for optically thick absorbers in the vicinity of 
1.8 < z < 4.0 foreground quasars. We presented 27 new 
quasar-LLS pairs with column densities 17.2 < log A^hi < 
20.9 and transverse distances of 22 kpc < R < 
1.7 h^^ Mpc from the quasars, corresponding to factors 
of guv = 5 — 8000 larger ionizing fiuxes than the am- 
bient UV background, if the quasars emit isotropically. 
The observed probability of intercepting an absorber is 
very high for small separations: six out of eight projected 
sightlines with transverse separations R < 150 kpc 
have an absorber coincident with the foreground quasar, 
of which four have A'hi > lO^^cm"^. The covering fac- 
tor of A^^Hi > lO^^cm-^ absorbers is thus ~ 50% (4/8) 
on these small scales, whereas < 2% would have been 
expected at random. 

Techniques for estimating quasar systemic redshifts 
from noisy rest frame UV spectra were presented and 
used to compute the systemic redshifts of the foreground 
quasars in our sample. An important area for future 
work is to obtain near-IR (7J-band) spectra of the z >2 
foreground quasars to determine more accurate redshifts 
from the narrow (cr < 200 km s"^) forbidden emission 
fines [O II] A3727 or [O III] A5007, which are the best 
predictors of systemic redshift. 

Much larger samples of optically thick absorbers near 
quasars are easily within reach, and can be compiled in a 
modest amount of o bserving on a 6- 10m class telescopes. 
iRichards et al.l l|200 41 constructed a photometric catalog 
of 2 < 3 quasars from the SDSS imaging data which 
has high completeness and efficiency, and photometric 
redshifts accurate to |Az| < 0.3. Indeed, many of the 
projected quasar pair s shown in Figure ^ w ere selected 
from this catalog (see iHennawi et aLll2006afl . The num- 
ber density of Zphot > 1.8 quasars in the photometric cat- 
alog is 10.5 deg^, so that the total number of projected 



sightfines expected with AO < 25" (i? < 150 h^^ kpc) 
for the current 8000 deg^ of SDSS imaging is 140, 
of which 70 would have super-LLSs if we extrapo- 
late from the ~ 50% covering factor of such systems 
observed in this study. At moderate resolution, suffi- 
cient SNR can be obtained in a 30 minute exposure at 
Keck, and about 2.5 times longer at the MMT. Thus, a 
7 night program on a 10m class telescope, or 18 nights 
at 6m class, targeting the closest projected pairs would 
result in a factor of about 20 more quasar-absorber pairs, 
with R < 150 kpc and logA'^ni ^ 19, than the four 
published here. If a multi-slit setup were used, wider 
separation projected pairs come for 'free', since these 
quasars can be targeted and observed simultaneously 
on the same masks, as was done for our Keck obser- 
vations. Finally, we note that a similar survey of pro- 
jected quasar pairs using the techniques presented here 
can be used to search for Mg II, C IV or other absorp- 
tion lines sytem near foreground quasar s (Bowei Tet alJ 
I2006bt iProchter. He nnawi. fc Prochas ka 20061 . 

A sample of ~ 100 LLSs near quasars would be of 
tremendous cosmological interest, providing new oppor- 
tunities to characterize the environments, emission ge- 
ometry, and radiative histories of quasars, as well as new 
laboratories for studying fiuorescent emission from op- 
tically thick absorbers and the physical nature of LLSs 
and DLAs. 
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Fig. 2. — Spectra of projected quasar pairs with optically thick absorbers coincident with the foreground quasar redshift. The left panels 
show the background (upper) and foreground (lower) quasars. Some of the metal line transitions discussed in Sl^and indicated by labels 
and dotted lines on the background quasar spectrum. Quasar emission lines are indicated on the foreground quasar spectrum. The panels 
on the right show a closeup of the H I Lya absorption in the background quasar spectrum. The (red) central curves indicate the best Voigt 
profile fit, (green) upper and lower curves show the fits corresponding to our ztlcr column density errors, and (blue) dashed lines indicate 
the redshifts of absorption components. The upper (purple) dotted line indicates the continuum level used for the fit and the lower (orange) 
dotted line indicates the noise level. The arrow indicates the foreground quasar redshift and the (gray) rectangle indicates our estimate 
of the redshift error. Three quasars are shown on the left for the triple systems SDSSJ0127+1507 and SDSSJ2347+1501. Cross hatched 
regions in the spectra of SDSSJ1204+0221 and SDSSJ2347+1501 indicate gaps in the spectral coverage of the Gemini GMOS spectrograph. 



QUASARS PROBING QUASARS 



15 




4000 



4500 



X (A) 



5000 



5500 



4180 



4200 




4000 



4500 



5000 



5500 



X (A) 



E 

05 




z, =3.39 

J0256+0039 

R =960 kpc/N 



3960 3980 4000 4020 4040 4060 
I0r 




4000 



4500 



5000 



5500 

X (A) 



6000 



6500 



7000 



5300 5320 5340 5360 5380 



Fig. 2. — continued. 



16 HENNAWI et al. 




4000 4500 5000 5500 6000 4500 4520 4540 

X (A) 




3600 3800 4000 4200 3620 3640 

X (A) 



Fig. 2. — continued. 



800 
600 
400 
200 

i2 600 
c 

o 

O 400 

200 




QUASARS PKOBiJNG QUAbAKS 




zz cn o m 03 
2 = = ^ ^ 


u 


1 INTO l\3 ^ _!. _!. - 






1 1 i-r- 

1 , : 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 

Ly« z,=2.18 
J081 4+3250 

R =61 kpc/h^jrf»M«*t'¥^^ 



17 



400 



300 



200 



100 




3800 



4000 



4200 

X (A) 



4400 



4600 



3860 



3880 





X (A) 




H 1 1 1 1 1- 



z, =3.20 

J0852+2637 



R =931 kpc/h 



CIV 



4220 4240 4260 4280 4300 
10r 




4000 4500 5000 5500 

X (A) 



6000 



6500 



5080 5100 5120 5140 5160 



Fig. 2. — continued. 



18 



HENNAWI et al. 




4000 



4500 



5000 



5500 

X (A) 



6000 



6500 




5260 5280 5300 



4000 




J I . I r, I I I I I L_ 



4500 



5000 



5500 3960 3980 4000 4020 4040 



X (A) 



E 

05 




4400 4600 

X (A) 



5200 4000 4020 4040 4060 
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APPENDIX 
UV ENHANCEMENT 
In this appendix we describe how we computed the quantity 
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Fig. 2. — continued. 



which is the maximum enhancement of the quasars' ionizing photon flux over that of the extragalactic ionizing back- 
ground, at the location of the background quasar sightline, assuming that the quasar emits isotropically. 



The ionizing photon number flux i^uvB (photons s ^ cm ^) due to the UV background is 



-FuvB 



-r-dv, 



(A2) 



where t'LL = 13-6 eV/ft. is the Lyman limit frequency, h is Planck's constant, and is the mean specific intensity of 
the UV background (erg s^^ cm^^ ster~^ Hz~^). For the intensity of the UV background, we use the version computed 
by F. Haardt & P. Madau (2006, in preparation), which considers the emission from observed quasars and galaxies 
after it is filtered through the IGM to yield the UVB as a function of redshift^^. 

The ionizing photon number flux from a quasar of speciflc luminosity (erg s~^ cm^^ Hz~^) at a proper distance 
rq is 

We assume that the quasars' spectral energy distribution obeys the power law form = L^^^^{v / vuS)'"^^ , blueward 



of i^LL- iTelfer et al.l l)2002|) measured an average power law slope of aq ~ 1.57 blueward of wavelengths < 1200 A, 
with the transition to this slope occurring in the range A — 1200 — 1300 A. 
The luminosity density at the Lyman limit for a quasar at redshift z is. 



(1 + ^) 



(A4) 



where c^l is the luminosity distance and f^^ is the specific flux of the quasar as observed on Earth at the redshifted 
frequency i/q = vu^/ {\ + z). We define the apparent magnitude at the Lyman limit, measured by an observer on Earth 

This model uses recent results for the quasar luminosity function, assumes that /esc = 0.1 of ionizing photons escape from galaxies, 
and assumes a power law index of a = 1.8 t o the s pectral energy distribution at wavelengths shortward of 912 AThe spectrum is available 
at http:/ / pitto.mib.infn.it / ~haardt / refmo del.html| 
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TABLE B2 

Tenative Optically Thick Absorbers Near Quasars 
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at the redshifted frequency according to 

/^^ = 10-°-4(48-60+™«i2) erg cni-2 Hz-^. (A5) 

For the majority of the quasars in Table ^ the Lyman hmit is below the atmospheric cutoff. Furthermore, even 
for higher redshifts, determining TOgi2 from our spectroscopic observations or the SDSS broad band magnitudes is 
complicated by the absorption from the Lyman-a forest and LLSs. Our goal is then to com pute m3i2 given th e SDSS 
apparent magnitude of a quasar, measured in the rest frame near-UV. To this end, we tie the lTelfer et alJ lj2002f) power 
law spectrum to the composite quasar spectrum of iVanden Berk et al.l 1)200 1|) at the wavelength 1285A, which is a 
clean region free of emission lines and in the wavelength range (A = 1200 - 1300 A) where l|Telfer et al.l|2003 ) observed 
the transition to a spectral index of ag. We then compute the 'k-correction' bet ween 771312 and the bl uest SDSS filter, 
which is also redward of the Lya line in the observed frame, by convolving the IVanden Berk et all l|2b01ll template 
with the SDSS filter curve. We work with a filter redward of Lya to avoid the Lya forest of the quasar, which is not 
correctly included in the quasar template. 

TENTATIVE OPTICALLY THICK ABSORBERS NEAR QUASARS 

Relevant quantities for a sample of tentative quasar-absorber pairs, for which we could not be sure that log A^hi > 17.2 
are given in Table lB2l Higher SNR and higher resolution spectra are required to make definitive conclusions about these 
systems, making them an interesti ng s et of targets for future research. Coordinates and SDSS five band photometry 
of these objects are given in Table [C4l of Appendix C. 

COORDINATES AND PHOTOMETRY 

Coordinates and SDSS five band photometry of the quasars in Tables ^ a-nd IB2I are listed below. We also list 
coordin ates and photometry for all five members of the projected group of quasars associated with SDSSJ 0127-1-1507 
in Table El 
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TABLE C3 

Coordinates and Photometry for Projected Quasar Pairs in Table 



Name 


RA (2000) 
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Note. — Coordinates, redshifts, and extinction corrected SDSS five band PSF photometry 
are given in the columns it, g, r, i, and z for the projected quasar pairs in Table The 
background and foreground quasars are labeled 'A' and 'B', respectively. 
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TABLE C4 

Coordinates and Photometry for Projected Quasar Pairs in Table IB21 



Name 


RA (2000) 


Dec (2000) 
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Note. — Coordinates, redshifts, and extinction corrected SDSS five band PSF p hoto metry 
are given in the columns u, g, r, i, and z for the projected quasar pairs in Table IB2I The 
background and foreground quasars are labeled 'A' and 'B', respectively. 



TABLE C5 

Coordinates and Photometry for, the Five Quasars In the Pro.iected Group 

SDSSJ0127+1507 



Name RA (2000) Dec (2000) 



SDSSJ0127+1508A 01:27:44.85 +15:08:58.0 2.60 20.95 20.56 20.50 20.38 20.33 

SDSSJ0127+1507B 01:27:42.57 +15:07:38.4 2.38 21.21 20.61 20.55 20.71 20.70 

SDSSJ0127+1506C 01:27:41.20 +15:06:55.9 2.18 18.73 18.51 18.36 18.36 18.17 

SDSSJ0127+1504D 01:27:40.77 +15:04:10.3 2.08 19.93 20.07 20.03 19.86 19.57 

SDSSJ0127+1506E 01:27:43.53 +15:06:48.4 1.82 21.21 21.25 21.14 20.94 20.99 



Note. — Coordinates, redshifts, and extinction corrected SDSS five band PSF photometry 
are given in the columns u, g, r, i, and z for the projected quasar pairs in Table The 
background and foreground quasars are labeled 'A' and 'B', respectively. 
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